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ABSTRACT 

We have performed a survey for substellar companions to young brown dwarfs in the Taurus star- 
forming region using the Wide Field Planetary Camera 2 on board the Hubble Space Telescope. In 
these data, we have discovered a candidate companion at a projected separation of 0'.'105 from one 
of the brown dwarfs, corresponding to 15 AU at the distance of Taurus. To determine if this object 
is a companion, we have obtained images of the pair at a second epoch with the adaptive optics 
system at Gemini Observatory. The astrometry from the Hubble and Gemini data indicates that the 
two objects share similar proper motions and thus are likely companions. We estimate a mass of 
5-10 Mj up for the secondary based on a comparison of its bolometric luminosity to the predictions of 
theoretical evolutionary models. This object demonstrates that planetary-mass companions to brown 
dwarfs can form on a timescale of t < 1 Myr. Companion formation on such a rapid timescale 
is more likely to occur via gravitational instability in a disk or fragmentation of a cloud core than 
through core accretion. The Gemini images also reveal a possible substellar companion (p — 0'.'23) to 
a young low-mass star that is 12'.' 4 from the brown dwarf targeted by Hubble. If these four objects 
comprise a quadruple system, then its hierarchical configuration would suggest that the fragmentation 
of molecular cloud cores can produce companions below 10 Mj up . 

Subject headings: protoplanetary disks — stars: formation — brown dwarfs — binaries: visual — 
stars: pre-main sequence 



1. INTRODUCTION 

Most of the known planets outside of the solar sys- 
tem have been detected through radial velocity measure- 
ments , and thus reside close to their par ent stars (a < 
5 AU. lMarcv et al.l2005HUdrv et al.l200l . Meanwhile, a 
small number of companions with masses below 10 Mj up 
have been detected at large separ ations through high- 
resolution imagin g (a ~ 30-300 AU iChauvin et al.l2 004: 
Lafreniere eFaLl 120081: iKalas et al.l 120081: iMarois et al.l 
2008) . The available data tend to favor core accretion as 
the mechanism responsible for the close-in giant planets 
(jSantos et al.|[2004TlFischer fc Valentil 2005). but gravita- 
tional instability in disks and fragmentation of molecu- 
lar cloud cores are plausible (perhaps likely) alternatives 
for producing the wide comp anions (|Lodato et alj|2005t 
iDodson-Robinson et al.ll2009f) . The formation models for 
wide planetary-mass companions are particularly in need 
of observational tests. 

Companions are expected to form much more quickly 
through disk instability and cloud core fragment ation 
than core accretion (r ~ 0.1 Myr vs. r ~ 5-10 Mvr. | Boss 
1998 1 [200l iBurkert fc Bodenheimer|[T99l iPoIiack et al. 
19961) . Thus, measurements of the frequency of plane- 
tary companions at very young ages can help identify 

1 Based on observations performed with the NASA/ESA Hub- 
ble Space Telescope and Gemini Observatory. The Hubble obser- 
vations are associated with proposal ID 11203 and were obtained 
at the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS 5-26555. 

2 Department of Astronomy and Astrophysics, The Penn- 
sylvania State University, University Park, PA 16802; 
kotl04@astro.psu.edu. 

3 Center for Exoplanets and Habitable Worlds, The Pennsyl- 
vania State University, University Park, PA 16802. 

4 Whitin Observatory, Wellesley College, Wellesley, MA. 



the dominant mode for the formation of these objects. 
To provide data of this kind, we have performed a search 
for substellar companions to 32 young brown dwarfs in 
the Taurus star-forming region (r ~ 1 Myr, d ~ 140 pc) 
using high-resolution images obtained with the Hubble 
Space Telescope. In this Letter, we present the discovery 
of a companion that has a mass of 5-10 Mj up and a pro- 
jected separation of 15 AU from one of the brown dwarfs 
in this survey. 

2. BACKGROUND 

The brown dwarf for which we have detected a sub- 
stellar companion is 2MASS J04414489+2301513 (hence- 
forth 2M J044144). It is a member of the Tau- 
rus st ar-forming regio n and it has a spectral type of 
M8.5 (jLuhmanl I2006D . which correspo nds to a mass 
of ~ 20 Mj UD for an ag e of 1 Myr ([Chabrier et al.l 
120001 ILuhman et all 120031 ). 2M J044144 exhibits mid- 
IR excess emission and strong Ha emission, indicat- 
ing the presence of a circumstellar d i sk and active 
accre tio n, respectively (ILuhmanl 120061 : ILuhman et al.l 
[20101) . IKraus fc Hillenbrand! (|2007D identified a pos- 
sible wide stellar companion at a separation of 12'.'4 
(a = 1700 AU at 140 pc). This candidate companion, 
2MASS J04414565+2301580 (henceforth 2M J044145), 
has been spectroscopically confirmed as a young low- 
mass star, and hence a likely member of Taurus (M 3- 
M4.5, ILuhman et aHl200ll IKraus fc Hillenbrand! l2009h . 

3. OBSERVATIONS 

3.1. WFPC2 Images 

We selected 2M J044144 as one of the targets for a sur- 
vey of late-type members of Taurus with the Wide Field 
Planetary Camera (WFPC2) on board the Hubble Space 
Telescope. The targets were imaged through the F791W 
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and F850LP filters with the PC camera within WFPC2, 
which contained a 800 x 800 CCD (0'.'046 pixel" 1 ). In 
each filter, we obtained two images at each position in a 
two-point diagonal dither pattern. For 2M J044144, the 
exposure times of the individual images were 260 and 
160 s for F791W and F850LP, respectively. The obser- 
vations of 2M J044144 were performed during one orbit 
on 2008 August 20. 

We red uced the WFPC2 images using the MultiDrizzlc 
software (Kockcm oer et all 120021) . To facilitate the de- 
tection of low-mass companions, we applied point spread 
function (PSF) subtraction to each br own dwarf in 
our su rvey using techniques described by iLuhman et all 
(2005). This process revealed a possible companion to 
2M J044144, which we refer to as 2M J044144 B. The 
images of this pair in F850LP before and after PSF sub- 
traction are shown in Figure [TJ We fit PSFs to both the 
primary and candidate companion to measure their flux 
ratios and astrometric positions in the two filters. In ad- 
dition to PSF subtraction of the primaries, we identified 
all point sources appearing in the WFPC2 images and 
measured aperture photometry for them. In the case of 
2M J044144, the PSF of the candidate companion was 
subtracted prior to the photometric measurement of the 
primary. The photometry for companion was computed 
by combining the aperture photometry of the primary 
with the flux ratios produced by the PSF fits. The errors 
in the photometry for the primary are dominated by the 
uncertainty in the correction for charge transfer efficiency 
while the errors in the measurements for the companion 
have a comparable contribution from the process of PSF 
fitting. The astrometric offsets and photometry for the 
two objects are presented in Tables Q] and [2] respectively. 
Additional details concerning the observations and data 
analysis will be provided in a forthcoming paper that 
presents the entire survey. 

3.2. NIRI+ALTAIR Images 

To constrain the relative proper motions of 2M J044144 
and its candidate companion, we o btained images of the 
pair at a second epoch using NIRI (IHodapp et alJ l2003) 
in conjunction with the ALTAIR adaptive optics system 
at the Gemini North telescope. NIRI was configured 
with the f/32 camera, which provided a plate scale of 
0'.'0214 pixel" 1 and a field of view of 22" x 22". Because 
natural guide stars were unavailable, ALTAIR was oper- 
ated in the laser guide star mode. The tip/tilt correction 
was performed with 2M J044145, which is the possible 
wide companion to 2M J044144 that was described in 
§ H We imaged 2M J044144 through the H and K' fil- 
ters using individual exposure times of 1 and 20 s. For 
each exposure time and filter, one image was taken at 
each position in a 3 x 3 dither pattern. After one set 
of dithered images was collected for all of the combina- 
tions of exposure time and filter, the entire sequence of 
observations was repeated. To provide a PSF for fitting 
the primary and candidate companion, we also observed 
2MASS 04383110+2310107. It was selected as the PSF 
star because it was similar to the science target in its 
sky position and the optical brightness, separation, and 
position angle of its tip/tilt star. 

The individual images were divided by flat field expo- 
sures taken with the Gemini Facility Calibration Unit. 
The resulting images in a given dither sequence were 



registered and combined, producing two final images of 
2M J044144 for each combination of exposure time (1 and 
20 s) and filter (H and K'). All of the combined images of 
2M J044144 and the PSF star exhibited FWHM~ 0"09, 
which is similar to the image quality from WFPC2. To 
measure the relative positions and flux ratios of the com- 
ponents of 2M J044144, we performed PSF fitting with 
the four long exposures in the same manner as done for 
the WFPC2 data. In Figure [T] we present one of the 
two long if'-band exposures before and after PSF sub- 
traction alongside the WFPC2 images. We also show 
a wider area of this NIRI image that encompasses the 
tip/tilt star, 2M J044145, which is resolved into a pair 
of objects that are separated by 0V23. We applied PSF 
fitting to the components of this pair as well to measure 
their flux ratios. This was done with the short exposures 
since the primary is saturated in the 20 s data. 

Our measurements of the astrometric offsets between 
2M J044144 A and B are listed in Table ffl We have 
adopted the average values produced by PSF fitting of 
the four long exposures. We have derived photome- 
try for the components of 2M J044144 and 2M J044145 
by combining the flux ratios from the PSF fitting with 
the unresolved H and K s data from the Point Source 
Catalog of the Two -Micron All-Sky Survey (2MASS, 
ISkrutskie et al.l [20061 ). We have assumed that the flux 
ratios at K s are the equal to the ratios that we have 
measured at K' . 

4. ANALYSIS 
4.1. Evidence of Binarity 

To investigate whether 2M J044144 B is a companion 
to 2M J044144 A, we begin by plotting a color-magnitude 
diagram in Figure [2] from photometry at F791W and 
F850LP for all unsaturated sources in our entire WFPC2 
survey of 32 brown dwarfs in Taurus. Four of the substel- 
lar primaries are absent from the diagram because they 
are saturated in one of the filters. The faintest primary is 
2MASS J04373705+233 1080, which is the co olest known 
member of Taurus (LP. ILuhman et al~ll2009f ). As shown 
in Figure [2j 2M J044144 B appears near the bottom of 
the sequence of the known members of Taurus, indicat- 
ing that it has the appropriate magnitudes and color for 
a low-mass companion. Given the spectral types of the 
known members, these data suggest a spectral type of 
M9.5-L0 for 2M J044144 B. In addition, the line connect- 
ing the primary and secondary is parallel to the cluster 
sequence, which is consistent with the coevality that is 
expected for the components of a binary system. 

We can use the data in Figure [2] to estimate the prob- 
ability that our companion survey would detect a field 
star that has the photometric properties of a low-mass 
Taurus member. We assume that the sequence of Taurus 
members becomes vertical at the faintest magnitudes in 
Figure [3] (771791 > 22), which is supported by the loca- 
tion of the coolest known member and the fact that the 
I — Z colors of field dwarfs do not increase significantly 
from M8 to mid-L ([Steele k, Howellsl [20001 iDahn et all 
I2002t iDobbie et al.ll2002[ ). In that case. ~ 7 objects in 
Figure [5] (in addition to the known brown dwarfs) have 
colors and magnitudes that are consistent with member- 
ship in Taurus. Therefore, the surface density of field 
stars that could be mistaken for Taurus members is ~ 7 



Planetary-mass Companion in Taurus 



3 



divided by the total area of the WFPC2 images, which is 
160 arcmin 2 . The probability of a field star of this kind 
appearing within O'.'l of any of the 32 primaries in our 
survey is ~ 10~ 5 . 

The two epochs of astrometry from WFPC2 and NIPJ 
provide an additional constraint on the companionship of 
2M J044144 B. We first consider whether the astrometry 
for 2M J044144 B is consistent with the negligible mo- 
tion expected for the vast majority of background stars. 
Accurate measurements of the proper motion for the pri- 
mary are not available. Therefore, we adopt for the pri- 
mary the average proper motion of the nearest group 
of Taurus members (/i a ,^s = +6.7, —17.7 mas yr _1 , 
iLuhman et al.l 120091) . Combining this proper motion 
with the astrometry from the first epoch, we find that the 
separation and position angle of the candidate compan- 
ion would have become 0'.'089 and 111.5° in the second 
epoch if it had no proper motion, which is inconsistent 
with our measurements in Table [1] 

According to the data in Table [lj 2M J044144 A 
and B maintained the same relative positions to within 
~ 10 mas over a period of 1.15 years, indicating that 
their proper motions agree at a level of ~ 9 mas yr _1 
in each direction. Only a small fraction of field stars are 
expected to have proper motions that are similar to that 
of 2M J044144 A. For instance, using multi-epoch Hub- 
ble images of Chamaeleon I that are comparable in depth 
to our WFPC2 data in Taurus, we find that ~ 0.3% of 
field stars have proper motions that are consistent with 
those of cluster members at the level of the uncertain- 
ties in the relative proper motions of 2M J044144 A and 
B. If we adopt this value for Taurus, which is similar to 
Chamaeleon I in its distance and the magnitude of its 
proper motion, then the probability of a field star ex- 
hibiting WFPC2 photometry that is indicative of a Tau- 
rus member, a separation of p < O'.'l from one of our 32 
primaries, and a proper motion similar to that of Taurus 
is - 3 x 10~ 8 . Therefore, we conclude that 2M J044144 B 
is a companion to 2M J044144 A. 

The possible wide companion 2M J044145 appears 
within one of the WFC arrays, but it is too heavily satu- 
rated for the detection of its candidate companion. As a 
result, we cannot place constraints on the proper motions 
of those two objects relative to each other. 

4.2. Physical Properties of Companion 

We can estimate the bolometric luminosity of 
2M J044144 B if we apply a bolometric correction to 
its photometry at K s . In previous studies of young 
brown dwarfs, we have adopted bolometric corrections 
derived for field dwarfs, but doing so may not be per- 
fectly valid for 2M J044144 B since some of the IR colors 
of young objects later than M9 differ fr om those of dwarfs 
(jCruz et alj|2009t ILuhman et alll2009D . To determine an 
appropriate value of BC#- for 2M J044144 B, we have es- 
timated bolometric luminosities for young brown dwarfs 
from M9.5-L0 for which relatively complete spectral en- 
ergy distributions have been measu red, consisting of the 
Taurus member KPNO 4 (M9.5, IBriceno et all [2001 
and the young field dwarfs 2MASS J0 1415823-4633574 
and 2MASS J02411151 -0326587 (L0, iKirkpatrick eTaTl 
120061: iCruz et all 12009). We have excluded from this 
sample brown dwarfs that exhibit disk emission in their 
mid-IR data or that have extinctions of Ay > 1. For 



each object, we constructed a spectral energy distribu- 
tion using a 0.8-2.5 /zm spectrum from the NASA In- 
frared Telescope Facility that was flux calibrated with 
2MASS photometry, 3.2- 9.2 /xm photometry from the 
Spitzer Space Telescope (|Luhman et all [2009. 20l3), a 
linear interpolation of the fluxes between 2.5 (im and 
3.2 fim, and a Rayleigh- Jeans distribution longward of 
9.2 /an. By integrating the flux in each distribution 
and combining the resulting luminosity with the K s 
photometry, we arrive at BCif = 3.42, 3.39, and 3.41 
for KPNO 4, 2MASS J01415823-4633574, and 2MASS 
J02411151— 0326587, respectively. These bolometric cor- 
rections are larger than those of dwarfs by ~ 0.2 mag 
(|Golimowski et al1l2004ft . Based on these measurements, 
we adopt BC^ = 3.4 for 2M J044144 B. 

The J- H color f r om 2M ASS and the 0.8-1.8 /jm spec- 
trum from|Luhman| ([20061 ) for 2M J044144 A+B indicate 
Ay ~ when they are compared to avera ge intrinsic col- 
ors a nd spectra of young brown dwarfs (jLuhman et all 
12010(1 . In the K band, the spectroscopy and 2MASS 
photometry of the system and our resolved photom- 
etry for the primary are slightly redder than a stel- 
lar photosphere with the spectral type of the primary, 
which is probably caused by emission from a circumstel- 
lar disk base d on the excess emiss ion detected at longer 
wavelengths ([Luhman et al.l 120101) . Therefore, we have 
not applied an extinction correction to the photome- 
try for 2M J044144 B. When we combine K s = 14.94, 
BC/f = 3.4, Mboi© = 4.75, and the average distance of 
Taur us (d = 140 pc, IWichmann et allll998[|Torres et all 
f2009l) . we derive log L hoi = -3.14 for 2M J044144 B. 

In Figure G2 we estimate a mass from the luminos- 
ity of 2M J044144 B by comparing it to the values pre- 
dicted by theoretical evolutionary models for an age of 
1 Myr, which is the average age for mem bers of Taurus 
and the isochronal age of the primary (jLuhman et all 
2009). Given the uncertainties in its age (few Myr) and 
luminosity (±0.12 dex), 2M J044144 B could have a 
mass between 5 and 10 Mj u p according to the models 
of lChabrier etaLl ([200(1 and Burrows et all (fl997l ). For 
comparison, we include in Figure [3] the luminosities of 
other resolved companions that appear to have masses 
below ~ 15 Mj up and that are younger than 10 Myr , 
which cons ist of 2M1207-3 932 B (IChauvin et all 12004ft. 
DHTau B (lltoh et alj2"005ft. CHXR 73 B (jLuhman et all 
l2006h. CT Cha B (ISchmidt et al.l l2008l). and IRXS1609 B 
(jLafreniere et al1l2008jl . 

5. DISCUSSION 

The first planetary-mass companion to be de- 
tected through direct imaging was 2MASS J12073346- 

3932539 B (a = 40 AU, M ~ 4 M 1lin 

I Chauvin et al.l l200l jBiller fc Close! [20071 : iGizis et all 
l2007tlDucourant et al.l l2008). which is bound to a young 
brown dwarf i n the TW Hya association (r ~ 8 Myr, 
M ~ 25 M 1 „ r . lGizisll2002ft . The new companion that we 
have presented, 2M J044144 B, is a younger analog of 
2MASS J12073346-3932539 B, thus demonstrating that 
wide planetary-mass companions to brown dwarfs can 
form on a timescale of r < 1 Myr. This result is consis- 
tent with the rapid formation expected from both grav- 
itational instability in disks and fragmentation of cloud 
cores. The latter is more likely to have produced these 
two secondaries given their relatively large masses com- 
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pared to the primaries (jLodato et al.l 120051) . Indeed, if 
2M J044144 A and B are members of a quadruple sys- 
tem, then its hierarchical configuration further suggests 
that 2M J044144 B formed by cloud core fragmentation. 
Additional data are needed to determine more defini- 
tively whether 2M J044144 A/B and 2M J044145 A/B 
comprise a quadruple system. 
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TABLE 1 

Astrometry for 2MASS J04414489+2301513 B 



Instrument 


9 


P.A. 


Date 




(arcscc) 


(dcg) 




WFPC2 


0.105±0.004 


120.4±2.2 


2008 Aug 20 


NIRI+ALTAIR 


0.108±0.005 


121.0±2.6 


2009 Oct 13 



TABLE 2 

Photometry for 2 MASS J04414489+2301513 and 2MASS J04414565+2301580 



2MASS 



01791 



nisso 



H 



J04414489+2301513 A 
J04414489+2301513 B 
J04414565+2301580 A 
J04414565+2301580 B 



18.93 ± 0.05 
21.16 ± 0.10 



17.85 ± 0.05 
19.91 ± 0.10 



13.94 ± 0.03 
15.62 ± 0.10 
10.17 ± 0.02 
13.03 ± 0.06 



13.40 ± 0.03 
14.94 ± 0.10 
9.94 ± 0.02 
12.59 ± 0.06 
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Fig. 1. — Left and middle: WFPC2 and NIRI+ALTAIR images of the young brown dwarf 2M J044144 before and after PSF subtraction 
(0'.'6 X (y/6). The candidate companion has the same position relative to the brown dwarf in each of these epochs, indicating that they 
have similar proper motions (Table[T]l- The separation of the companion is 0"105, which corresponds to 15 AU at the distance of Taurus. 
Right: An expanded version of the NIRI+ALTAIR image (13" X 13") showing both 2M J044144 and the young low-mass star 2M J044145. 
A candidate substellar companion to the latter is detected in these data (p = ty.'23). The components of 2M J044144 and 2M J044145 may 
comprise a young, low-mass quadruple system. 



Planetary-mass Companion in Taurus 



16 



24 



• 


I I 


*■ .' .« 

#■* % •• 


• • 

• • 

A ® 


• . .V. • '.. * 

■ * — *V I % ■ . 

r >o ^ »". J- * . 

< v.f.-j'^A • • 
^.r/.W.'.'f '.:*« • •. . 


• 

• • 

B 

• 






, ■, , r- ■, ; 





0.5 1 1.5 



Fig. 2. — Color-magnitude diagram constructed from WFPC2 images of 32 young brown dwarfs in the Taurus star-forming region. We 
show the targeted brown dwarfs that have unsaturated photometry (large points) and all other point sources in these images (small points). 
The components of 2M J044144 are indicated (open circles, Figure [TJ. The secondary appears near the bottom of the sequence of known 
substellar members of Taurus, which is consistent with a young low-mass companion. 
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Fig. 3. — Luminosity estimates of 2M J044144 B and other young low-mass companions (M < 15 Mj up ) compared to the luminosities 
as a function of age predicted by the theoretical evolutionary models of Chabricr ct al. (2000, dashed lines) and Burrows ct al. (1997, 
dotted lines). CHXR 73 B and CT Cha B have nearly identical positions in this diagram. These models indicate a mass of 5—10 Mj up for 
2M J044144 B. The luminosities of all of these companions have uncertainties near ±0.12 dex. 



